Measuring entanglement in condensed matter systems 
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We show how entanglement may be quantified in spin and cold atom many-body systems using standard 
experimental techniques only. The scheme requires no assumptions on the state in the laboratory and a lower 
bound to the entanglement can be read off directly from the scattering cross section of Neutrons deflected 
from solid state samples or the time-of-flight distribution of cold atoms in optical lattices, respectively. This 
removes a major obstacle which so far has prevented the direct and quantitative experimental study of genuine 
quantum correlations in many-body systems: The need for a full characterization of the state to quantify the 
entanglement contained in it. Instead, the scheme presented here relies solely on global measurements that are 
routinely performed and is versatile enough to accommodate systems and measurements different from the ones 
we exemplify in this work. 



I. INTRODUCTION 



Interacting quantum many-body systems generally exhibit 
correlations between its constituents. At sufficiently low tem- 
peratures, near the ground-state, these correlations possess 
quantum mechanical features, namely entanglement. Com- 
pared to its classical counterpart, entanglement is extremely 
complex. Its full characterization generally requires the mea- 
surement of a number of observables that grows exponentially 
with the number of constituents of the quantum many-body 
systems. On the one hand, the ability to create entangle- 
ment merely by cooling an interacting quantum many-body 
system provides the attractive opportunity of using this en- 
tanglement to carry out quantum information processing tasks 
such as measurement based quantum computation [ 1 1 or adi- 
abatic quantum computation |2| which gain their power ex- 
actly because of the complex structure of entanglement. How- 
ever, the very same setting offers significant challenges, as it 
is much harder to analyze theoretically and, crucially, experi- 
mentally: While for the few-particle systems that can now be 
prepared in highly controlled environments such as ion traps, 
it is possible to fully characterize the state in the laboratory 
by quantum state tomography [3|, the situation in condensed 
matter systems is far more challenging. Firstly, the number 
of subsystems tends to be much larger, the level of control 
over states and Hamiltonians is more restricted and crucially 
the available measurements for condensed matter systems are 
much less general: Local measurements addressing individ- 
ual constituents are usually not available and one has to rely 
on global measurements such as those obtained in scattering 
experiments to draw conclusions about the system. These are 
of course by no means sufficient to fully characterize the state 
in the laboratory. How might one still be able to say some- 
thing about the entanglement that is available? One approach 
could be to, e.g., model the system with a certain quantum 
Hamiltonian and compare it to a classical model. If the pre- 
dictions from the quantum model matches the measurement 
results while the classical does not, and the simulated quan- 
tum state displays entanglement, one might conclude that the 
state in the laboratory is indeed entangled. This however is 
a fallacy. Consider the following example of two spins (4): 
Suppose one measures the correlation (a\<J 7 ^) — (af) (erf) and 



obtains the result —1. This measurement is consistent with 
both the maximally entangled state \ip) = (|t|) - ||t))/\/2 
and the separable state q = (|ti)(ti| + l4-t>Utl)/ 2 - Hence, 
without further assumptions, one may not decide whether the 
state in the laboratory is entangled or not. A possible assump- 
tion may be that the system is in thermal equilibrium at some 
known temperature and that the Hamiltonian that governs the 
system is known precisely. But obtaining knowledge about 
the Hamiltonian experimentally is even harder than to obtain 
the state itself! A technique to decide without a doubt wether 
entanglement is contained in a given system should hence not 
rely on knowing the Hamiltonian, it should, in fact, not rely on 
any kind of knowledge about the system other than measured 
data, but be able to quantify entanglement by just relying on 
measured observables. 

Here we present a scheme to quantify entanglement in con- 
densed matter systems that fulfils all the above requirements 
and relies only on measurements that are already available: 
Neutron scattering from spin systems and time-of-flight imag- 
ing of cold atoms. Hence, we show that it is possible to 
directly — without any assumptions — measure entanglement 
in many-body systems. To this end we exploit the substantial 
body of work concerning the characterization (which states 
are entangled), quantification (how much of it do we have) and 
verification (on the basis of simple measurements we need to 
answer the previous questions) of entanglement (see for a 
tutorial review and [6 1 for a advanced and very comprehensive 
review) that has been established in quantum information sci- 
ence. More precisely, we combine methods for determining 
the presence of entanglement in quantum many-body system 
121 with proposals for the quantitative verification of entan- 
glement in few body quantum systems J4] l9 fflD to achieve 
an experimentally accessible method for measuring entangle- 
ment. 



II. SPIN SYSTEMS: NEUTRON SCATTERING 

One of the standard tools to analyze condensed matter sam- 
ples is neutron scattering, see, e.g., Ref. fl2l . The deflected 
neutrons carry information about both the structural and mag- 
netic properties of the sample, which can be read of the dif- 
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ferential scattering cross section lfl3l . The Fourier transform 
of the magnetic cross section gives access to spin correlations 
in reciprocal space such as the positive hermitian observable 

S(q) = Y, a =x, y ,z S<*(q), where 



S a (q) 



1 



M ^ 

»,3 



(1) 



= ( ? o )>*i = ( ? ~o )• = ( o -i ) are Pauli spin ma- 
trices acting on lattice site i located at r^, and M is the total 
number of spins, so the number of lattice sites. Analogous 
observables can be obtained for spin systems realized in ion 
traps or cold atoms in optical lattices, in which one can di- 
rectly measure spin correlations by light scattering [ 14]. 

In the following we will show how a measurement of ob- 
servables of the type in Eq. ([TJ alone is sufficient to quantify 
the entanglement contained in the sample. We will consider 
systems comprised of spin 1/2 particles (generalizations to 
higher spins are entirely straightforward) on a lattice and set 
out to derive a lower bound to the entanglement that is consis- 
tent with the measurement data. It will turn out that it is possi- 
ble to find a lower bound that is a simple function of the static 
structure factor (S(q)). Hence, for this, no assumptions about 
the state are necessary, in particular, no knowledge about the 
Hamiltonian is required. 

In the following we present definitions and a derivation that 
will lead to the central result in Eq. (|6j, which provides a lower 
bound on the entanglement of any spin-state under investiga- 
tion that can be used directly on experimental data. Several 
entanglement measures may be expressed in the form 131 IT31 



£c(q) = max'j 0, 



mm 

wewnc 



(2) 



where W denotes the set of hermitian operators that fulfil 
(W) > for separable states (i.e., the set of entanglement 
witnesses |[T6lD and C distinguishes the quantities: E.g., if C is 
the set of operators W fulfilling (W) < 1 for separable states 
then £c is the robustness of entanglement 1 17 1, for C = {W £ 
W 1 1 — W > 0}, £c measures the generalized robustness of 
entanglement (TBI, and if C = {W £ W 1 1 + W > 0} then 
£c is equal to the best separable approximation |[T9l . In fact, 
Ee n , m , C„, m := {W e W | - nl < W < ml}, is an en- 
tanglement monotone for every n, m > |fT5l . One can now 
exploit the fact that for any choice W <E W H C one obtains a 
lower bound to £ c , 

£c{0) > max|o, -tr[#p]| for all W € W n C. (3) 

Given this expression, it is possible to arrive at lower bounds 
to £q by simply constructing operators W £ W n C that are 
functions of observables that are within experimental reach. 
This works of course for any spin system and any observable. 

In the following, we will focus on the best separable ap- 
proximation and the observable S. Consider the operator 
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W(q) - -S(q) - 1, 



(4) 



E{q) 
1- 



0.8 
0.6 
0.4- 
0.2- 




r/2 



r/4 



t/ 2 3tt/4 



1y 



FIG. 1 : Lower bound E on the entanglement — as measured in terms 
of the best separable approximation (BSA) — for a thermal state q = 
exp(-/3H)/Z, /3 J = J/k B T = 1, of the Heisenberg model in 
Eq. (jTJ. For every q, E(q) provides a lower bound to the BSA. 
The square lattice with open boundary conditions and lattice constant 
a — 1 has 30 x 30 lattice sites and (S(q)) was obtained using the 
generalized directed loop quantum Monte Carlo algorithm |22| of 
the ALPS package (23). Inset shows E( , =| ), n = 0, 4, 6, 7 (top 
to bottom), as a function of the temperature. Lines are a guide to the 
eye. Note that the BSA is upper bounded by unity, a bound that E(q) 
saturates at low temperatures and q = 0. 



for which we now show that W € WnC. We find W(q)+1 > 
and for a product state q = ®iQi, we have 



2M 



(5) 



which is non-negative as the last term is and the first term 
may be bounded by using the uncertainty relation 5^„(1 — 
(a?) 2 ) > 2. Hence, W E WnC, i.e., for every state q and 
every q, the quantity 

E(q) = max jo, 



1 - 



1 

2M 



(6) 



provides a lower bound to the A/ -partite entanglement (as 
measured in terms of the best separable approximation) con- 
tained in q. Similar bounds may be derived for all entangle- 
ment measures that fall into the general framework of Eq. Q. 

At this point, we would like to emphasise again that E(q) 
gives a lower bound to the entanglement for any state on the 
lattice - irrespective of how it has been prepared, what the 
temperature is, or what the underlying Hamiltonian of the sys- 
tem might be. 

As an example, we consider thermal states of the antiferro- 
magnetic Heisenberg model 



a=x,y,z 



(7) 



on a square lattice. This model has been analyzed in great 
detail in the literature (see Ref. [20] for a review) using sev- 
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eral analytical and numerical techniques. In Ref. Il2p the 
two-dimensional Heisenberg antiferromagnet copper deutero- 
formate tetradeurate has been analyzed experimentally using 
extensive neutron scattering measurements and it has been 
suggested — under the assumption that the system is indeed 
described by the Heisenberg Hamiltonian with known cou- 
pling constants and by a comparison of a classical to a quan- 
tum description — that entanglement is present in this system. 
Using E(q), the presence of entanglement can not only be 
confirmed but, as it is a lower bound to the best separable ap- 
proximation, also quantified. In Fig.[T] we show E(q) for a 
thermal state of the Heisenberg model at different values of 
J/3 as obtained from a quantum Monte Carlo simulation us- 
ing the generalized loop algorithm [22 1 of the ALPS library 
ll23l . The plot shows that entanglement is present up to fairly 
high temperatures, i.e., measuring entanglement is well within 
experimental reach (in GTl . e.g., the sample was at a temper- 
ature of 1.5 K and the data well fitted by J = 6.19 meV, i.e., 
ksT/J = 0.02). In addition, the plot exemplifies the qual- 
ity of our bound: The best separable approximation is upper 
bounded by unity and at low temperatures E(0) s» 1. E(q) 
also scales properly with the system size: For the ground state 
of the Heisenberg model it is known that S(q) ~ \q\ for small 
|q| E0), i.e., E(q) - 1 ~ \q\, see also footnote EU. 



where Uj(x,y) = Jdz/^ff), f iti (x,y) =: f(x,y) . We 
now set out to show that 



E(x,y)=m3 X {o,{N)-j^\} 



f(x,y) 



(11) 



provides a lower bound to the entanglement contained in the 
state in the laboratory, which constitutes the main result of this 
section. Here, (N) — I s me expected total number 

of atoms. 

As we are concerned with massive particles, we will, in the 
following restrict the state space to states g that have a finite 
mean number of particles, tr[giV] < 00, and commute with the 
particle number operator N. In other words, we are concerned 
with states respecting the particle-number superselection rule 
(SSR) - the only physical states allowed in this setting of in- 
distinguishable massive particles [30). These states are of the 

form g = Y,n=o p nqPn =■ 0w=o &v> where P N projects 
on the sector with constant particle number N. The SSR also 
restricts the allowed physical operations to operations com- 
muting with N [ 3 1 1 . Consider now 
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00 
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tr[&vW]}, 



(12) 



III. BOSONS IN OPTICAL LATTICES: TIME-OF-FLIGHT 
DENSITIES 

A standard measurement in the context of ultracold atoms is 
the following: One switches off all potentials, allows the atom 
cloud to expand freely and then takes an absorption image 
of the atoms, which reveals the velocity, quasi-momentum, or 
time- of -flight distribution of the atoms before the expansion. 
This technique was used to show Bose-Einstein condensation 
into the lowest-momentum state (25], to demonstrate the Mott 
insulator - superfluid transition of bosons in optical lattices 
ll26l . and to observe Fermi surfaces of fermions in optical lat- 
tices ll27l to name just a few. We focus on the situation in 
which bosons of mass m are kept in an optical lattice with 
lattice constant a. After a time of flight t, the density of the 
atoms reads (see, e.g., Refs. [28, 29]) 



n(r) 



i,3 



fi-j 



(8) 



where 



f iJ (k) = (f t r\w(k)^e 



w(k) is the Fourier- transform of the Wannier function cen- 
tred at zero, and bi annihilates a boson at site i located at r*j. 
The resulting absorption image is then the integral along the 
optical axis, say, the z-direction, of this density, i.e., 

n{x,y) = ^fc,j{x,y)(b\b i ) =: (h(x,y)), (10) 

i,3 



where Cjv is the set of hermitian operators W acting on the 
subspace of constant particle number N that fulfil cN ± W > 
for some constant c > independent of N and tr[gArV[^] > 
for separable qm- In the Appendix we show that £ is an en- 
tanglement monotone under LOCC operations that preserve 
the total number of particles, i.e., that commute with iV (and 
hence it is also an entanglement monotone under SSR-LOCC 
operations - LOCC operations that preserve the local parti- 
cle number; for a discussion of entanglement under SSR see 
Ref. 1301 "). 

We now show that W N = P N (n/f(x 1 y) - N)P N e C N . 
To this end let be a state vector on the subspace of constant 
particle number N. Then, with M being the number of lattices 
sites, we find 



l i3 



(13) 



'■3 

^E 

1,3 



NM, 



i.e., MN ± Wn > 0. Furthermore, for separable qn, we find 
for i ^ j that tr[^StS 3 -] = and hence 

ti[g N n] = f(x, y) ^ tr^St^] = f(x, y)N, ^ 



i.e., tr^TvVKjv] > 0. Hence, E C^r, which implies that 

- min \x[q n W] > \i[g N (N - n/f(x, y))] (15) 

W£C N 
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FIG. 2: Lower bound E(x, y) on the entanglement for a thermal 
state g = exp(—H /ksT)/Z with constant filling factor (hi) — 1 
of the Bose-Hubbard model in Eq. l |16[ >. The three-dimensional cubic 
lattice with periodic boundary conditions and lattice constant a — 1 
has 10 x 10 x 10 lattice sites and {n(x, y)) was obtained using the 
same numerical code as for Fig. [T] Left plot shows E(x, y) as in 
Eq. {HI for PU = 1/5, J/U = 0.01. Right plot shows ff), 
n = 64, 48, 44, 36, 34, 33 (top to bottom) as a function of the tem- 
perature (black) and of the tunnelling amplitude J/U (gray). Lines 
are guides to the eye. 



and thus, for all x, y, the quantity E(x, y) in Eq. (Hi provides 



a lower bound to the M -partite entanglement (as measured in 
terms of £) available in the system. Here, n(x, y) = (h(x, y)) 
is obtained in standard time-of-flights measurements and E is 
a lower bound for any state on the lattice. 

As an example, we consider thermal states of the Bose- 
Hubbard model on a three dimensional cubic lattice, 



H 



U 



^2n,i(hi - 1) - ju^', 



(16) 



where hi — b\bi, (•, •) denotes summation over nearest neigh- 
bours, J accounts for tunnelling of atoms between adjacent 
sites, U is the strength of the on-site repulsion of atoms, and 
the chemical potential /i controls the particle number. Bosons 
in deep optical lattices are (up to an harmonic confinement) 
well described by this model, which displays a quantum phase 
transition from a Mott insulator (small J/U) to superfluid 
(large J/U) that was observed in [26 1 via the interference 
pattern displayed in (h(x,y)). In Fig. [2] we show E(x,y) 
for a thermal state of the Bose-Hubbard model as obtained 
from a quantum Monte Carlo computation using the ALPS li- 
brary 1231 . We can see that E(x,y) increases linearly with 
J/U and stays finite up to high temperatures. Furthermore, 
in experiments the atoms are harmonically trapped leading to 
an extension of the atoms over the lattice that is on the order 
of the system size considered here. Hence, quantifying entan- 
glement in these systems is already well within experimental 
reach. 



IV. CONCLUSION 

We have derived lower bounds to the entanglement con- 
tained in lattice systems. These lower bounds are a simple 
function of routinely measured observables and do not require 
any additional information about the system. This makes not 



only the verification but also the quantification of entangle- 
ment in condensed matter samples possible. In other words, 
without making any assumptions (such as the temperature, the 
Hamiltonian that governs the system, the way the state was 
created), entanglement can be directly measured using only 
measurements that already belong to the toolbox for the anal- 
ysis of quantum many-body systems. The presented schemes 
straightforwardly generalize to other many-body systems and 
observables and we foster the hope that they will see direct 
application to experimentally realized situations and inspire 
further generalizations. 
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VI. APPENDIX 

We set out to show that on the set of states g that respect 
the particle number superselection rule (i.e., are of the form 
8 = E??=o PnqPn) and have a finite mean number of par- 
ticles, tr[giV] < oo, the quantity £ as in Eq. 



12 1 is an entan- 



glement monotone under LOCC operations that preserve the 
total particle number. First, we note that £ is well defined: The 
set Cm is a subset of hermitian operators on A/jv, which we de- 
note by H n, and the trace is understood to be over A/jy. Then, 
for given N and qn the mapping Hn ^'[qn 1 ] is contin- 
uous, i.e., the minimum is attained on Cn as Cm is compact: 
Cm = An n Sn, where 



A 



N 



w en 



N 



cN <W <cN\ 



(17) 



is compact and the complement of 



S N = \w eH N H[qnW] > for separable^} (18) 

is open in Hn- Let W € Hn but not in Sn. Then there exists 
a non-separable qn such that tr[^ArVl^]/2 = — e < 0. Then 
for all W' G H N with \\W' - W\\ < e, we have tr[g N W'} = 
tr^ivC^' - W)\ + tx[g N W] < \\W' - W\\ - 2e < -e < 0. 
Hence, the minimum is attained and, for given qn denoting 
the minimizer by Wg N , we may write 

E(g) = max jo,- tx[g N W SN }\. (19) 
^ jv=o ' 



Finally, the series is absolutely convergent: 

- cNti[g N ] < tr[g N W] < cNtr[g N ] 



(20) 
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for all W G C N , i.e., |tr[gjvV^]| < cNtr[g N ] for all W G Cat, 
and therefore 

oo oo 

53 MetfW 6lf ]\ < c Nti[g N ] = ctx[gN] < oo. (21) 



7V=0 



N=0 



Hence, £ is well defined. 

Now let A k , k = 1, 2, . . . , be Kraus operators of the form 
A k = 2 jv=o -^n (i- e -' we onr y a ll° w operations that preserve 
the total number of particles) with ^ fc AtAfc < 1 and A k — 
® s A k s , where the direct product refers to some partition of 
the system. For a given state g = (J)^ =0 £>n> denote 



gk 



Pk 



N=0 



N=0 



Then, with K' = {k G K | J2n=o *[(&Wp } < 0}, 



(22) 



R k = A k gA\ = A k N g N {A%)\ 

N=0 

Pk — trfiife], If = {k \pk > 0}, and for k E K write 

oo oo 

RklPk = - AW(iw) f == Qn- (23) 



which is finite. Hence, we may interchange the sums to find 



keK 



PkE{g k ) 



OO 

E 

JV=0 



tr[^(^(i^,)t^i 



(27) 



< - 53 tr [^^«] < ^(e), 

where (*) holds if JJfcGiC' (^iv) ^jv * s an e l ement of 

Cn, which we now set out to show. As J^ k A\A k < 1, we 

also have Ylki^kV — ^> i- e -' bounding positive semi- 
definite operators by zero, 



cNi ± 53 [AWwp Ab 



keK' 



> J2(A k N )^cNA k N ± 53 (A%^w 6% Al 

k keK' 

> 53(i^) t (ciVl±W^)i^>0, 



(28) 



keK' 



keK 



keK 



53 PkE{Sk) = 53 Pfe max \ °' - 53 tr $V^e 

oo 

53 Pk 53 tr^Wj* ] 

oo 

53 53 tr[^(i^)t#^ife 



keK' N=0 



(24) 



We already know that X^jvLo Mf^v^s* ] converges abso- 
lutely. We also have 

53 Pfcltrl&l^j, ]| - 53 HA k N Q N (A k N )^W §% )\ 

keK' keK' 

<cN^[A k N gN(A k N )^} (25) 
fe 

= ciVtr[^(53(i^)ti^)], 
which is upper bounded by cJVtr[gjv], i.e., 

oo oo 

53 53 PkH9%W d% }\ < c 53 Na[g N ] = ctr[gN], (26) 



N=0keK' 



N=0 



i.e., J2 k eK'(A%yWgk A k N G An- Now let gN be separable. 



Then 



tr 



QN 



(53(i^)t%i^ 



feeK' 



53 tr[i^(i^)t^ 



(29) 



where each summand is non-negative: The A\\gN{A%y are 
separable as the A k N are local operations and gN is separable. 



Therefore, as Wzk G Sn we have 



A k N g N (A k N ^W A 



> 0, 



(30) 



i.e., finally, EkeK'i^VW^M G 5 
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